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Abstract. We study theoretically the measurement of a mechanical oscillator using a 
single two level system as a detector. In a recent experiment, we used a single electronic 
spin associated with a nitrogen vacancy center in diamond to probe the thermal motion 
of a magnetized cantilever at room temperature {Kolkowitz et al., Science 335, 1603 
(2012)}. Here, we present a detailed analysis of the sensitivity limits of this technique, 
as well as the possibility to measure the zero point motion of the oscillator. Further, 
we discuss the issue of measurement backaction in sequential measurements and find 
that although backaction heating can occur, it does not prohibit the detection of zero 
point motion. Throughout the paper we focus on the experimental implementation of 
a nitrogen vacancy center coupled to a magnetic cantilever; however, our results are 
applicable to a wide class of spin-oscillator systems. Implications for preparation of 
nonclassical states of a mechanical oscillator are also discussed. 
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1. Introduction 

Recent interest in mechanical oscillators coupled to quantum systems is motivated 
by quantum device applications and by the goal of observing quantum behavior of 
macroscopic mechanical objects. The past decade has seen rapid progress studying 
mechanical oscillators coupled to quantum two-level systems such as superconducting 
qubits [U [2j [3], and single electronic spins [4], and theoretical work has explored strong 
mechanical coupling to collective atomic spins [6]. Recently, it was proposed that 
a mechanical oscillator could be strongly coupled to an individual spin qubit [7J [8]. 
Experiments based on single spins coupled to mechanical systems have demonstrated 
scanning magnetometry [9], mechanical spin control [10], and detection of mechanical 
motion [HI [12]. In parallel, pulsed spin control techniques have attracted renewed 
interest for decoupling a spin from low frequency noise in its environment, extending 
its coherence [13], while also enhancing the sensitivity of the spin for magnetometry 

[HE5J[T6J[I7]. 

In this paper we consider pulsed single spin measurements applied to the detection 
of mechanical motion at the single phonon level. We extend the analysis presented in 
our recent work [12] , providing a detailed theoretical framework and a discussion of 
measurement backaction. The central concept of our measurement approach is to apply 
a sequence of control pulses to the spin, synchronizing its dynamics with the period of 
a magnetized cantilever, thereby enhancing its sensitivity to the motion. By measuring 
the variance of the accumulated phase imprinted on the spin by the oscillator during a 
measurement, we directly probe the average phonon number, despite the fact that the 
oscillator position is linearly coupled to the transition frequency of the spin. We derive 
the conditions for observing a single phonon using the spin as a detector, and find that 
these conditions coincide with that of large effective cooperativity, sufficient to perform a 
two-spin gate mediated by mechanical motion [18]. Further, we consider the backaction 
arising from sequential measurements and show that this does not prohibit single phonon 
resolution. Throughout the paper, we focus on the specific spin-oscillator system of a 
magnetized cantilever coupled to the electronic spin associated with a nitrogen- vacancy 
(NV) center in diamond. For realistic experimental parameters we find that this system 
can reach the regime of large cooperative spin-phonon coupling, and the spin may be 
used to measure and manipulate mechanical motion at the quantum level. 

We begin in Sec. |2] by introducing the coupled system and spin control sequences, 
and calculate the signal due to thermal and driven motion of the oscillator. Then in 
Sec. [3] we derive the optimal phonon number sensitivity, and show the relation between 
strong cooperativity and single phonon resolution. Finally, in Sec. H] we consider the 
limit of zero temperature and calculate the signal due to zero point motion, including a 
discussion of backaction heating for sequential measurements. 
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2. Coherent sensing of mechanical motion 

2.1. Model 

We consider the setup shown schematically in Fig. (TJ in which a magnetized cantilever is 
coupled to the electronic spin of a single NV center. The magnetic tip generates a field 
gradient at the location of the NV, and as a result its motion modulates the magnetic 
field seen by the spin causing Zeeman shifts of its precession frequency. To lowest order 
in small cantilever motion, the precession frequency depends linearly on the position of 
the tip and is described by the Hamiltonian (h = 1) 

H = -a z + ^ (a + a 1 ") a z + H OSCJ (1) 

where a z is the Pauli operator of the spin and a is the annihilation operator of the 
oscillator. For a spin associated with an NV center in diamond, we take |t) = \m s = 1} 
and \\) = \m s = 0) in the spin-1 ground state of the NV center, and safely ignore the 
\m s = — 1) state assuming it is far detuned by an applied dc magnetic field. A is the 
detuning of the microwave pulses used for spin manipulation, which plays no role in what 
follows and we take A = throughout the paper. The spin-oscillator coupling strength 
is A = g e fiBG m xo/h, where g e ~ 2 is the Lande ^-factor, /ie is the Bohr magneton, 
G m is the magnetic field gradient along the NV axis, and xq = J%~/ 2muo is the zero 
point motion of the cantilever mode of mass m and frequency ujq (we included % in the 
definitions of A and xq for clarity). The damped, driven oscillator is described by 

H osc = u a j a + + H Ar , (2) 

where 

#7 = 9k (a + a f ) (b k + b\)+J2 Ukb\b k (3) 

k k 

describes dissipative coupling to a bath of oscillators b^, characterized by damping rate 7 
and temperature T. Finally, ifdr describes a coherent oscillator drive which we consider 
briefly in Sec. 12.41 Note that in Eq. ([1]) we have temporarily omitted intrinsic spin 
decoherence due to the environment; we will include this explicitly in Sec. [3j 

2.2. Spin echo and multipulse sequences 

The motion of the oscillator imprints a phase on the spin as it evolves under Eq. ([T|), 
which can be detected using spin echo [31 H5] . or more generally a multiple pulse 
measurement. Throughout the paper we focus on Carr-Purcell-Meiboom-Gill (CPMG) 
type pulse sequences, consisting of equally spaced n pulses at intervals of time r, as 
depicted in Fig. [TJ After initialization in a it/ 2 pulse prepares the spin in an 
eigenstate of & x , \i[) ) = |(|t) + ID) with (t[) \ a x |^ ) — 1- The spin is then allowed to 
interact with the oscillator for time t, accumulating a phase, and during which time we 
apply a sequence of n pulses which effectively reverse the direction of spin precession. 
At the end of the sequence, a final n/2 pulse converts the accumulated phase into 
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Figure 1. (a) Schematic of the setup. A single spin can be used to measure mechanical 
motion via magnetic coupling, (b) Toggling sign of the interaction describing 7r pulses 
flipping the spin. Each sequence begins and ends with n/2 pulses, and 7r pulses flip the 
sign of the interaction at regular intervals of time r. Thin dashed line shows oscillator 
position, which is synchronized with pulse sequence for r = tt/luq as shown. The total 
sequence time is t = 2r for spin echo and t = Nt for CPMG. 



a population in |f) which is then read out. By applying both initial and final 7r/2 
rotations about the same axis, we measure the probability to find the spin in its initial 
state |^o) at the end of the sequence, given by 

P(t) = \(l + (a x (t))), (4) 

where angle brackets denote the average over spin and oscillator degrees of freedom. 
Our choice to measure a x probes the accumulated phase variance; this is crucial for our 
purpose because the average phase imprinted by an undriven fluctuating oscillator is 
zero. In contrast, by applying the first and final ir/2 pulses about orthogonal axes one 
would instead measure a y , which probes the average accumulated phase. 

The sensitivity of the spin to mechanical motion is determined by the impact of 
the oscillator on the spin coherence (a x (t)). The key to maximizing this impact is to 
synchronize the spin evolution with the mechanical period using a CPMG sequence of 
7r pulses, increasing the accumulated phase variance and improving the sensitivity as 
discussed in the context of ac magnetometry [20]. Choosing r = tx/ujq between the 7r 
pulses, we flip the spin every half-period of the oscillator and maximize the accumulated 
phase variance. At the same time, these pulse sequences decouple the spin from low- 
frequency magnetic noise of the environment, extending the spin coherence time T 2 
[T5l IT6] . We describe the effects of the applied n pulses using a function f(t, r), which 
flips the sign of the spin-oscillator interaction at regular intervals of time r as illustrated 
in Fig. [TJ In this toggling frame, the interaction Hamiltonian is 

H int (t) = h z X(t)f(t,r), (5) 

where X = a + & and X(t) = e lHoBct Xe~ lHoBat . We calculate the spin coherence, 
(a x (t)) = (u^{t)a x U{t)), where the evolution operator is U(t) = Te~^ z So dt'x(t')f(f,T) 
and T denotes time ordering. Since the interaction is proportional to a z , it leads to 
pure dephasing and we obtain [2T] 

(a,(t)) = (re-^ 2 Te-^ 2 ) osc , (6) 
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where we used (a x (0)) = 1, the average (■) BC is over oscillator degrees of freedom, T 
denotes anti-time ordering, and the accumulated phase operator is 

= A f dt'X(t')f(t',r). (7) 
jo 

The spin coherence in Eq. ([6]) can be calculated using a cumulant expansion, which is 
vastly simplified by noting that the full Hamiltonian in Eq. (CD), including the oscillator 
drive and ohmic dissipation, is quadratic in X. As a result, the second cumulant — which 
in general corresponds to a Gaussian approximation — in the present case constitutes the 
exact result. We use this below to calculate the coherence for both thermal and driven 
motion. 

Another consequence of the fact that H is quadratic in X is that the effect of 
the pulse sequence is completely characterized by its associated filter function [201 122] , 



F(ut) = with f(u>) = J dte f(t,r). The filter function describes how two- 

time position correlations (x(t)X(t')^ of the oscillator affect the spin coherence in the 
second cumulant in the expansion of Eq. fl6]). For the pulse sequences illustrated in 
Fig. [H the corresponding filter functions are 

n/ \ f 8sin 4 (wr/2) spin echo 

F(ujt) = < „ (8) 

[ 2 sin 2 (AW/2) [1 - sec (cur/2)] 2 CPMG 

Note that phase-alternated versions of CPMG, such as XY4, which vary the axis of 7r 
pulse rotation in order to mitigate pulse errors, are also described by the above model 
in the limit of ideal pulses. 

2.3. Thermal motion 

As discussed above, the spin coherence in Eq. ([SD is given exactly by its second order 
cumulant expansion. Since the total sequence time is t — Nr, the coherence depends 
only on the time r between tt pulses, 

(a x (t = Nr)) = e-^ T \ (9) 

where 

Xn(t) = A 2 / ^ISxiu), 10 

and S x (u) = I dte iujt \({X{t), X(0)}) is the symmetrized noise spectrum of X. For the 
damped thermal oscillator described by H osc in the abscence of a drive, the symmetrized 
spectrum is (ks = 1) 

= , , 2w 7wcoth(u;/2T) 
(to 2 — coq) 2 + 7 2 ur 

where 7 = uq/Q is the mechanical damping rate due to coupling to the ohmic 
environment at temperature T. 

We plot the spin coherence due to thermal motion in the classical limit T ^> u>o in 
Fig. [21 The impact of the oscillator is greatest when the pulse sequence is synchronized 
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Figure 2. Spin coherence for CPMG sequence with N = 8, with undriven thermal 
oscillator at temperature T = 10 oj (a) and T = 1000 w (b) and values of Q shown. 
Solid lines show full spin coherence with collapses and revivals, and dashed lines show 
oscillator- induced dcphasing resulting in envelope decay given by Eq. (|12j) . Here we 
took X/loq = 0.01 and neglected intrinsic spin decoherence, T\ = T 2 — > oo. 



with the cantilever frequency, r = (2k + 1)ix/ujq with k an integer. At times r = 2/c7r/u;o, 
the accumulated phase due the oscillator cancels within each free precession time, so that 
the accumulated phase variance averages nearly to zero and the coherence revives. We 
stress that this structure of collapse and revival can arise from purely classical motion; 
it is simply a consequence of averaging the phase variance accumulated by the spin 
over Gaussian distributed magnetic field fluctuations with a characteristic frequency. In 
addition to collapses and revivals, the finite Q of the cantilever also causes dephasing 
of the spin which leads to an exponential decay factor of the envelope as e _r ^ T '. In the 
limit Q > 1 and T > ojq, the dephasing rate is given by 

r> ~ 3iVy 7 frith + I) , (12) 



2, 

where r\ = \/ujq is the dimensionless coupling strength and n t h = (e w °/ T — l)" 1 is 
the thermal occupation number of the oscillator. We provide a derivation of Eq. ( IT2|) in 



Appendix A| Increasing Q not only increases the depth of the collapses in spin coherence 
due to the oscillator, but also decreases the overall spin dephasing resulting in more 
complete revivals, as shown in Fig. [2j We also see that increasing the temperature 
increases both the depth of collapse and the dephasing. Below in Sec. [3] we use these 
results to calculate the lowest temperature motion that can be detected, characterized 
by the phonon number sensitivity at the optimal pulse timing r = ti/uq. 

2.4- Driven motion 

It is straightforward to include the effects of a classical drive through H^ r in Eq. (JT]). 
This simply adds a classical deterministic contribution to X(t), and we can decompose 
the accumulated phase in Eq. (0) it as <f> = 0dr + 0th where 

dr = XAj dt cos (u t + 6 )f(t, t) (13) 
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is the classical accumulated phase due to the drive. Here, A is the dimensionless 
amplitude of driven motion and 8 is its phase at the start of a particular measurement. 
We assume that the cantilever drive is not phase-locked to the pulse sequence, so 6q is 
random and uniformly distributed between and 2tt. Using Eq. fl6]) and averaging over 
60 we obtain 

(a x {t = NT)) = J [a{T)]e-x»W (14) 

where J is the zeroth order Bessel function [15], a(r) = t]A^J2F(uot), and Xn(t) is 
the thermal contribution given by Eq. fTiTJ]) . For a strong drive, thermal fluctuations 
are unimportant and the signal is given by the Bessel function. For a weak drive, 
comparable to thermal motion with \A\ 2 ~ n t h, both thermal and driven contributions 
may be important as illustrated in Fig. |3]and observed in experiment [12]. In Fig. [3] we 
see that, unlike thermal motion (see Fig. [2]), driven motion can lead to dips in the spin 
coherence below zero. In the remainder of the paper we focus on detecting undriven 
thermal or quantum motion with the drive switched off. 




Figure 3. Spin coherence from combined thermal and driven motion for drive 
amplitudes shown. For a weak drive, both driven and thermal contributions are 
important. The dips in the spin coherence below zero arise from driven motion, 
described by the Bessel function in Eq. (Tl4| . Parameters are u>o/2tt = 1 MHz, 
T = 50 w , Q = 100. 



3. Phonon number sensitivity 

In this section we discuss the sensitivity limits of the spin used as a detector of undriven 
mechanical motion. By comparing the signal from thermal motion to the relevant noise 
sources, we obtain the phonon number sensitivity. We then discuss the sensitivity in 
several limits relevant to experiments. 

3.1. Signal 

The impact of an undriven thermal oscillator on the spin coherence in a spin echo 
or CPMG measurement sequence is described by Eqs. (jl]) and ([9]). In addition to 
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its coupling to the oscillator, the spin is also coupled to an environment which leads 
to intrinsic decoherence and degrades the signal. For an NV center, decoherence or T 2 
processes are caused by a 1 % natural abundance of 13 C nuclear spins in the otherwise 12 C 
lattice. Flip-flop processes between pairs of these nuclear spins produce low frequency 
magnetic noise which leads to decoherence of the form e~ N( - T ^ T ' 2 ^ for a CPMG sequence 
with N pulses [201 122] • Note that T 2 here refers to the decoherence time in a spin 
echo sequence (i.e. N=l), typically ~ 100 fis in natural diamond and up to ~ 2 ms 
in isotopically pure diamond [23]. An added benefit of multipulse sequences is the 
enhanced spin coherence time, T 2 = A r2//3 T 2 , due to dynamical decoupling [13]. Finally, 
spin-lattice relaxation due to phonon processes leads to exponential decay on a timescale 
Ti, typically ~ 1 ms at room temperature and up to ~ 200 s at 10 K [21]. Including 
these intrinsic sources of spin decoherence, as well as the oscillator-induced decoherence 
given in Eq. ( Tl2|) . the probability to find the spin in its initial state given in Eq. (J4j) 
is modified as 

P(t = Nt) = \ (l + e-^+^^e-V) _ S(t), (15) 
where we have isolated the coherent signal due to the oscillator, 

S(r) = i e -^+C^) 3 ) e -IV (l _ e -(xN(r)-T,r)\ (16) 

Note that we have accounted for the oscillator-induced decoherence T^t which 
diminishes the coherent signal we are interested in. 

We can obtain a simple analytic expression for the signal in the limit Q ^> 1. In 
this limit the oscillator spectrum is well- approximated by Lorentzians at w — ±a>o, 

(w-w ) 2 + 7 2 /4 (w + w ) 2 + 7 2 / 4 
Using Eq. (fTTj) with Eq. ( fit)]) we obtain a compact analytic expression for \ with no 
further approximation, which we provide in Appendix A We choose the pulse timing 
t to maximize the impact of the oscillator motion on the spin coherence, providing 
optimal sensitivity. This is achieved by setting r = tt/ujq, flipping the spin every half 
period of the oscillator and resulting in the maximum accumulated phase variance. For 
N ^> 1, the filter function with r = tt/uq is well- approximated by a Lorentzian centered 
at ojq of bandwidth buo/N, where b ~ 1.27. Together with Eq. (1171) this yields 

WrfQN (_ 1\ 
XMvr/wo) ^ 7z ■ ,^ /An [rhh + 77 > ( 18 ) 



7T (1 + bQ/N) V 2 
and substituting this into Eq. f ll6p we obtain the signal. 



3.2. Sensitivity 

To find the sensitivity we must account for 
photon shot noise into a single parameter 
measurements is a = l/K^/M, where M = 



noise. We combine spin projection and 
K so that the noise averaged over M 
ttot/Nr is the number of measurements 
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of duration Nt that can be performed in a total time t tot . It follows that the minimum 
number of phonons that we can resolve in a given time t tot is 

° 

^min = I jo / j- f' I 

\dS/dn th \ 



and the corresponding phonon number sensitivity is £ = riminVnot- Using Eqs. ( 1161) and 
(j!8p with r = 7i /uq we obtain 

£ ~ ^ e"l"* (l + ^ (20) 
where we have expressed the total spin dephasing in terms of a single pulse number, 



71 / 7T \ 3 37T?7 2 /_ 1 



(21) 



which combines both intrinsic and oscillator-induced decoherence. Eqs. ( 1201) and (121j) 
reflect the competition between the oscillator damping rate 7 = ouq/Q, the intrinsic 
decoherence times T\ and T 2 of the spin, and the measurement bandwidth bu^/N. It is 
clear from Eq. ( fl8l) that increasing the number of pulses increases the coherent signal 
due to the oscillator; however, this also leads to increased spin decoherence. As a result, 
the resolvable phonon number is minimized at an optimal number of pulses, 



iVopt = N^-bQ + ^Nl + GbQN^ + (bQ) z . (22) 

Note that the optimal pulse number is always set by the spin decoherence, iY opt ~ N^, 
with only a prefactor of order one depending on Q. Neglecting pulse imperfections, the 
optimized sensitivity is determined by an interplay of Q, T\ and T 2 in Eq. (120|) . In 
practice, the optimal pulse number may be very large due to long spin coherence times, 
and pulse errors may play a role as discussed further below. 

In Fig. H]we plot the sensitivity as a function of pulse number N, and the optimized 
sensitivity as a function of coupling strength A. To check the validity of the above 
approximations it is straightforward to calculate the phonon number sensitivity directly 
from Eqs. ( fXOT) and (lllj) . The numerically exact sensitivity is shown in Fig. H] in 
agreement with our analytic results. In the remainder of this section we discuss the 
sensitivity in several experimentally relevant limits. 



3.3. Optimal sensitivity and cooperativity 

An important limit for current experiments is one where the spin coherence is much 
longer than the oscillator coherence during the measurement, corresponding to N<j, > Q. 
We assume that the spin coherence is dominated by intrinsic sources described by T\ and 
T 2 , and that the oscillator-induced spin decoherence Y ^ can be neglected, well-justified 
in the limit of weak coupling. Within these limits, the optimal number of pulses is 
Aopt ~ ^4> an d the optimized sensitivity is 

7T 3 / 2 
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(a) (b) 




N v = X/u 

Figure 4. (a) Phonon number sensitivity £ versus the number of pulses N for values 
of Q shown and A/wo = 0.01. Lines show the analytic result in Eq. ([20]) and points 
show the full numerical result using Eqs. (fTU| and (fTTj) . Squares mark the sensitivity 
at the optimal pulse number N opt . (b) Sensitivity optimized with respect to N, versus 
the coupling strength X/uig for the same values of Q as in (a). Squares mark the 
optimized sensitivity at X/loq = 0.01, corresponding to the squares in (a). The dotted 
lines mark a sensitivity of £ = l/VHz. Parameters in both plots are loq/2-k = 1 MHz, 
T 2 = 100 fjs, T x = 100 ms, T = 4K and K = 0.3. 



where the cooperativity is 

C = (24) 

7 

and T 2 = N^T 2 is the enhanced spin coherence time due to decoupling. For a large 
number of pulses, the enhanced spin coherence N 2 / 3 T 2 may be very long, and ultimately 
the spin coherence may be limited by T\ which is not suppressed by decoupling. In this 
case Eqs. (123]) and (j24j) are simply modified by T 2 — > T\. The cooperativity parameter C 
is ubiquitous in quantum optics, and marks the onset of Purcell enhancement in cavity 
quantum electrodynamics. In the present case, C > 1 is the requirement for a single 
phonon to strongly influence the spin coherence, leading to a measurable signal despite 
the relatively short coherence time of the oscillator. The condition C > 1 to resolve a 
single phonon can be simply understood: if the spin coherence is much longer than the 
oscillator coherence, i.e. Q N^, the accumulated phase variance increases at a rate 
~ A 2 /7 (see Eq. ( fl8l) with sequence time Nr ~ N/luq) and the maximum interrogation 
time (assuming that oscillator-induced decoherence is negligible) is T 2 . 

With feasible experimental parameters, T 2 ~ T\ ~ 10 ms, \/2n ~ 150 Hz, 
uq/2it ~ 1 MHz and Q ~ 1000, a cooperativity of C ~ 1 can be reached. In current 
experiments, NV centers exhibit a 30% contrast in spin-dependent fluorescence, and 
collection efficiencies of 5% are realistic [20l [25] . These parameters yield K ~ 0.3 and 
an optimal phonon number sensitivity of £ opt ~ l/v^Hz with N ~ ~ 15000 pulses. 
Due to long spin coherence times T\ and T 2 , the optimal pulse number may be very 
large, and in practice finite pulse errors may play an important role in limiting the spin 
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coherence. For example if the number of pulses is limited to N ~ 1000, a sensitivity of 
£ ~ 3/ \/Hz can be reached. We discuss this further below when we calculate the signal 
due to zero point motion. 

3.4- Ideal oscillators and ideal spin qubits 

While the cooperativity regime describes an important part of parameter space, it is 
useful to briefly consider two more simple limits that describe features in Fig. HJ First, 
we consider a harmonic oscillator that remains coherent for a much longer time than 
the entire pulse sequence, satisfying Q ^> N. In this limit, the long oscillator coherence 
time plays plays no role and the optimal sensitivity is limited only by the spin coherence, 
£o P t ~ 1 / (K A 2 T| ^Jujq/N). This limit can be seen on the left side of Fig. Hk, where the 
sensitivities for different values of Q fall on the same curve at low pulse numbers N. 

Finally, we consider the limit of very strong but incoherent coupling where the spin 
decoherence is dominated by the oscillator, i.e. T$ becomes larger than 1/T\ and I/T2. 
This limit is reached when either the intrinsic spin decoherence is negligible or for very 
strong coupling, rfn^ ^> Q/fwo^) 3 , Q/^qTi- In this limit, the coherent signal is large 
due to strong coupling, but saturates at a low number of pulses; further increasing the 
coupling strength only increases the oscillator-induced decoherence, reducing the signal. 
This is reflected in Fig. |4b, where we see that increasing the coupling strength larger 
than if > l/^n^T\ no longer improves the optimized sensitivity but instead degrades 
it. 

4. Detecting quantum motion 

Above we found that for realistic experimental parameters, a single phonon can be 
resolved in one second of averaging time. This raises the intriguing question of whether 
a single spin can be used to sense the quantum zero point motion of an oscillator in its 
ground state. It also implies that we must consider the effect of measurement backaction, 
which we have so far ignored in our discussion. To address these questions we analyze 
the experimentally relevant scenario where the spin is used to detect the motion of a 
mechanical resonator which is externally cooled close to its ground state. 

4-1. Measuring a cooled oscillator 

Even at cryogenic temperatures, a mechanical oscillator of frequency ojq/2-k ~ MHz has 
an equilibrium occupation number n t h much larger than one. For this reason we assume 
that the mechanical oscillator is cooled from its equilibrium occupation n t h to a much 
lower value Uq ~ 1 using either optical cooling techniques [26] or the driven spin itself 
[3 [27] . An important consequence of cooling below the environmental temperature is 
the effective reduction in Q of the oscillator. For an oscillator coupled to both a thermal 
environment and an external, effective zero temperature source for cooling, the mean 
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phonon number satisfies 




(25) 



(26) 



and in order to maintain n < 1 we require 7 coo i > 7n t h- As a result, the relevant 
decoherence rate of the oscillator is the rethermalization rate 771th • For this reason, 
to calculate the signal from a cooled oscillator we replace the equilibrium thermal 
occupation number nth by the effective occupation no — > in all expressions, while 
at the same time replacing the intrinsic Q by the reduced, effective quality factor 
Qcs = ^o/7cooi ~ Q/nth- 

4-2. Single shot readout 

In Sec. [3] we calculated the sensitivity £, which reflects the minimum detectable phonon 
number n m [ n that can be resolved in one second of averaging time. For the following 
discussion it is useful to convert the sensitivity to a minimum detectable phonon number 



per single measurement shot, n m i nj i = £ / \j (Nir / ujq) , where we have taken the total 
measurement time to be t tot = Nr and r = ti/uq. Assuming single shot spin readout 
(K — > 1), which has been demonstrated at low temperature [25], and using Eq. ( J20l) we 
obtain 



where C e s = A 2 T 2 /7n t h is the reduced, effective cooperativity. We see that under 
the assumption iV ~ A^ opt ^> Q c q, the ability to resolve ground state fluctuations of 
a cooled oscillator within a few spin measurements requires C c s > 1, which is the 
same strong cooperativity condition required to perform a quantum gate between two 
spins mediated by a mechanical oscillator [TS]. Alternatively, n min l corresponds to the 
occupation number required to produce a signal S of order one in Eq. (TIB"]) . It provides a 
convenient way to directly compare the sensitivity with the backaction due to sequential 
measurements, as discussed below. 

In Fig. [5] we plot the calculated signal due to zero point motion, assuming that 
the mechanical oscillator is cooled near its ground state n = and using the reduced 
quality factor Q e s- These plots show that the intrinsic coherence times typical for NV 
centers are more than sufficient to resolve single phonons provided enough pulses can be 
applied to exploit the full spin coherence. In practice, the limiting factor is likely to be 
finite pulse errors, which limit the absolute number of pulses that can be applied before 
losing the spin coherence. To estimate the effect of finite pulse errors, we include the 
calculated signal assuming additional spin decoherence of the form e~ N ^ Nc with a cutoff 
pulse number N c . Pulse numbers of N ~ 160 have been demonstrated in experiment 
[13] . and with further improvements this can be increased to more than N ~ 1000. 





(27) 
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Figure 5. (a) Spin coherence with uq ~ for increasing pulse number and Q c h = 100 
with A/cj = 0.01. (b) Optimal signal as defined in Eq. (|T6| from zero point motion. 
Solid lines show optimal signal assuming unlimited pulse number, while dashed lines 
include a simple treatment of pulse errors with N c = 1000 as described in the text. 
Parameters are T 2 = 100 [MS, T x = 100 ms, lo /2it = 1 MHz. 

Based on this we plot the modified signal using N c ~ 1000 and find that even with a 
limited number of pulses, zero point motion results in a significant signal for realistic 
coupling strengths. 

4-3. Backaction 

The result that a single spin magnetometer can resolve the quantum zero point motion 
of a mechanical oscillator calls for a discussion of measurement backaction. We begin by 
noting that, despite the linear coupling of the spin to the oscillator position in Eq. (TfJ, 
the described measurement protocol is sensitive to the variance of the accumulated phase 
~ (X 2 ), which we obtain by averaging independent spin measurements. As a result, 
our approach does not correspond to standard continuous position measurement [28J, 
nor does it implement a quantum nondemolition measurement of the phonon number, 
since the interaction in Eq. (JTJ) does not commute with n. In principle, by cooling 
between measurements our approach may be used to measure the phonon number with 
arbitrary precision. Nonetheless, the effect of the spin's backaction on the oscillator is 
both a practical issue and interesting in itself, and could be used to prepare nonclassical 
mechanical states. We describe two possible approaches to observe the influence of 
measurement backaction on the oscillator. 

First, we consider directly probing the projective nature of the measurement. 
For simplicity we assume that the oscillator is initially in its ground state and 
decoupled from the environment, and assume single shot spin readout. In a single 
measurement sequence, the oscillator experiences a spin- dependent force according to 
Eq. (j5J). Measuring (a x ) = ±1 at the end of the sequence projects the oscillator onto a 
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superposition of coherent states [HI 129] . 

. , . \ia) ± I— ia) , , 

17 V2 ± 2e~ 2a2 

where a = NX/2uq is the total displacement for a sequence of N ^> 1 pulses and 
r = tt/loq. The probabilities to measure |±) are given by 

P± = \(l±e-^), (29) 

which shows, consistent with the discussion above, that for a measurement strength 
a > 1 the oscillator in its ground state can significantly affect the spin dynamics. 
To observe the backaction of this measurement on the oscillator we can perform a 
second spin measurement, which is sensitive to the state of the oscillator conditioned on 
the first measurement. In principle, by using techniques developed in cavity quantum 
electrodynamics, this procedure can be used to fully reconstruct the conditionally 
prepared oscillator state [30] . 

Let us now consider an alternative, indirect way to observe backaction by 
performing many successive measurements. Again beginning with the oscillator near 
its ground state, the first measurement projects the oscillator into one of the states 
|"0±)- By averaging over the two possible spin measurement outcomes, the resulting 
mixed oscillator state is 

p osc =p + |<0+) (lp+\ +P- {lp-\ , (30) 

and we see that on average the oscillator energy has increased by \a\ 2 . Repeating 
this measurement many times, without cooling between measurements, the oscillator 
amplitude undergoes a random walk of stepsize ±a, and on average the phonon number 
increases approximately linearly in time. This corresponds to backaction heating 
described by an effective diffusion rate, 

Cb ,= ^. (3!) 

Combining the measurement backaction with intrinsic mechanical dissipation and 
external cooling, the average occupation number satisfies 

(h) = -(7 + 7cooi) (n) + 7f2 t h + -Dba, (32) 
and for 7 cool 3> 7 the steady state phonon number added due to backaction is 
£ba NQ eSV 2 

n ba = = — -. • (33) 

7cooi 4vr 

We see that increasing the coupling strength not only improves the single shot resolution 
^min,i) but also leads to backaction heating of the oscillator. For sufficiently strong 
coupling, the steady state backaction phonon number nba exceeds the phonon number 
resolution, and the inferred phonon number is determined by backaction. We thus take 
the sum n mcas = n min 1 + as a measure of the minimum inferred phonon number. 
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77 = X/iUQ 



Figure 6. Solid lines show total inferred phonon number given by Eq. (|33[) from 
combined phonon resolution and backaction heating. Dashed lines show sensitivity 
and heating contributions. For each value of Q e ff we set N = Q c ff/5. 



Note that for simplicity in this discussion we have assumed the limit N <C Q e s, in which 
the oscillator is coherent within each measurement sequence. Within this limit we find 

_ ira NQ cS 7] 2 
rVeas- 8 ^ 2iV2 + 4?r • (.M) 

The total inferred phonon number n meas is shown in Fig. |6] as a function of the coupling 
parameter rj and a fixed number of pulses N = Q c s/5. In this case n meas is minimized 
for rj ~ 1/VQcfr, where it reaches a value of 77. mcas ~ Observing this minimum 

in the phonon number resolution as a function of coupling strength would provide an 
indirect signature of measurement backaction. This observation may be more feasible in 
near term experiments than directly observing projective backaction as discussed above. 



5. Summary and conclusions 

We have presented the sensitivity limits of a novel position sensor consisting of a 
single spin. For realistic experimental parameters, we predict that a single NV 
center in diamond can be used to resolve single phonons in a cooled, magnetized 
mechanical cantilever. The condition to resolve single phonons is that of strong effective 
cooperativity, the same condition needed to perform a quantum gate between two spins 
mediated by a mechanical oscillator. For even stronger coupling, the backaction of 
the spin on the oscillator can be probed directly or indirectly, and used to prepare 
nonclassical mechanical states. 
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Appendix A. Analytic signal for thermal motion in high Q limit 

Here we sketch the derivation of Eqs. (1121) and (TT8]) . The impact of the oscillator on the 
spin coherence is given by Eqs. ffTOj) and ffTTj) . 

, 2 fduF(uT) wcoth(W2T) /A , 

XN (r) = auwA" / -^ w -^±-L- r (A.1) 

To perform this integral it is useful to decompose the filter function as 



N-l 



F(cor) = l-cos(AW)+ J2 C 1 ~ cos^s^-j (1 - cos^)) 

j=0 



(A.2) 



where Sj = (j + l/2)r and tj = (N — j)r. We first consider the high temperature limit, 
T ^> Uq, in which we can approximate coth(u;/2T) m 2T/u. The result is a sum of 
integrals of the form 

lrr , l9 f doj 1 — cos (cut) 9/ . . . , . „, 

4Tu; o7 A 2 / 1 ; 2 = 77 2 2n th g t , A.3 

which can be done exactly. In the limit Q 3> 1 we obtain 

g (t) = 7 t + (l - e -^/ 2 cos(w *)) - p-e~ yt/2 sm(co t). (A.4) 

To calculate we need the spin coherence at the revivals, given by Xn(j = ^n/uo). 
To first order in ft, we have q(t = Att/uq) ~ 3yt/2, and to this order the only nonzero 
term in Eq. ( 1A.1I) is due to the 1 — cos(Nojt) terms in Eq. (1A.2I) . The result is Eq. (I12j) . 



Next we derive Eq. f fT8|) at r = vr/wo in the limit N ^> 1. Here we use 
the fact that the filter function near r ~ tt/uq may be rewritten for > 1 as 
F(u7i/u ) ~ 2A r2 sinc 2 [7rA r (a; — w )/2], and in turn this function is well-approximated 
by its Lorentzian envelope, 

(w - uj y + {bu Q /Ny/A 
where we obtain the effective bandwidth buo/N by fitting the extrema of Eq. (1A.5|) to 
a Lorentzian which yields b ~ 1.27. At the collapse time r = n/uo, we can approximate 
$x(u>) by the Lorentzian spectrum given in Eq. (jT7|) . Using Eqs. (fTTj) and (]A.5j) . x(tt/u)q), 
the integrand is simply the product of two Lorentzians and performing the integration 
yields Eq. (jig). 
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